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We study the micromechanics of coUagen-I gel with the goal of bridging the gap between theory 
and experiment in the study of biopolymer networks. Three-dimensional images of fluorescently 
labeled collagen are obtained by confocal microscopy and the network geometry is extracted using 
a 3d network skeletonization algorithm. Each fiber is modeled as a worm-like-chain that resists 
stretching and bending, and each cross-link is modeled as torsional spring. The stress-strain curves of 
networks at three different densities are compared to rheology measurements. The model shows good 
agreement with experiment, confirming that strain stiffening of collagen can be explained entirely 
by geometric realignment of the network, as opposed to entropic stiffening of individual fibers. 
The model also suggests that at small strains, cross-link deformation is the main contributer to 
network stiffness whereas at large strains, fiber stretching dominates. Since this modeling effort uses 
networks with realistic geometries, this analysis can ultimately serve as a tool for understanding how 
the mechanics of fibers and cross-links at the microscopic level produce the macroscopic properties 
of the network. While the focus of this paper is on the mechanics of collagen, we demonstrate a 
framework that can be applied to many biopolymer networks. 



Collagen is the most abundant animal protein 1 and 
its mechanics have been studied in great detail 2^. It 
takes on many morphologies, including skin, tendons, lig- 
aments, individual fibers, and gels. Of particular interest 
is the mechanics of coUagen-I gels, shown in Figure [T^. 
These gels provide a relatively simple structure that can 
be noninvasively observed by confocal microscopy [5] and 
used as a scaffold for growing artificial tissues |4], and as 
a 3d environment for studying cell motility [5] and tumor 
invasion [Sl|7]. A critical first step in understanding these 
systems is to develop a model for the collagen gel alone. 
In this paper we give a successful theoretical model of 
the micromechanics of realistic networks. 

CoUagen-I gels belong to a class of materials known 
as biopolymers. Other examples include actin, found in 
the cytoskeleton and fibrin, a component of blood clots. 
A common feature of biopolymer networks is their abil- 
ity to strain stiffen by 2-3 orders of magnitude at large 
strains (Figure llj) . The cause of this strain stiffening is 
not well understood. Storm et al. ^ attributed strain 
stiffening in all biopolymer networks to the pulling out 
of entropic modes of individual filaments. Their cal- 
culation required the assumption that deformations are 
afFine. Later, Heussinger et al. [HI [10] showed how one 
could deconstruct the network deformation into a set of 
floppy modes. They concluded that accounting for the 
non-affinity was necessary in describing the elastic prop- 
erties of the network. Onck and colleagues [TTJ [12] have 
proposed the alternative hypothesis that strain stiffening 
is due to the rearrangement of the fibers as the network 



is strained (Figures [Tli and [T];) . Resolving this debate 
has been difficult since almost all theoretical analysis has 
been on artificially generated networks. The few exam- 
ples of quantitative comparisons to experiment in the lit- 
erature [HI [131 [H] are not able to quantitatively fit the 
full stress-strain response of the gel at varying densities 
using a single set of parameters. 

In this paper, we bridge the gap between model and 
experiment. Three dimensional images of fluorescently 
labeled collagen gels at different densities are imaged by 
confocal microscopy (Figure lltL) and the network geome- 
try is extracted using a custom FIbeR Extraction (FIRE) 
algorithm (Figure ^p) [3] . The gel is modeled as a ran- 
dom network of cross-linked fibers, as described below, 
and the stress-strain response is compared to that mea- 
sured by an AR-G2 rheometer. Good agreement between 
model and experiment is obtained by fitting a single pa- 
rameter, the cross-link stiffness. The experiments are 
described in detail in [3]. 

In the model for the collagen gel, each fiber is treated 
as a discrete worm-like-chain (WLC) [11] which resists 
stretching and bending, and each cross-link is treated as 
a torsional spring, thus more stiff than a freely rotat- 
ing pin joint but less stiff than a welded joint of fixed 
angle. The stretching modulus of an individual fiber is 
given by Kg — EA, where E is the Young's modulus and 
A is the cross-sectional area. The Young's modulus of 
a fiber in aqueous conditions has been estimated to be 
between 30-800 MPa [m [HI [18] and we use a modu- 
lus of 50 MPa, which fits the data well and is also close 
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FIG. 1: A typical gel (1.2 mg/ml, 25.6 ^m x 25.6 /im x 25.6 
/Ltrri). a) Maximal Intensity Projection along the axis that is 
perpendicular to the focal plane of the microscope, b) Projec- 
tion of 3d network extracted by FIRE, c) Reduced network, 
where elements that don't contribute to the network stiffness 
have been removed for improved computational efficiency, d) 
Deformation after 50% tension, e) Deformation after 50% 
shear, f) Comparison of the stress-strain response between 
the model and experiment. 



to the value chosen by Styhanopoulous et al. (79MPa) 
to fit their model [H]. It has been shown that a single 
fiber will stiffen by a factor of 2-4 when strained [TH] . We 
choose here to use a constant E both to reduce the num- 
ber of parameters in the model and to see if geometric 
reorientation of the network is enough to explain strain 
stiffening. Stylianopoulos and Barocas 14] also explored 
the bilinear and exponential constitutive relations for the 
individual fibers and observed only minor effects on the 
macroscopic network behavior. The radius of each fiber 
is r = 30 nm [3, 19 . The bending modulus of the fiber is 
given by Kb ^ EI ^ 32 pN-^m^, where / = 7rr'*/4 [20]. 
No cross- linking agent has been added to the gel and very 
little is known about the nature of the naturally formed 
collagen cross-links. We find that we can fit all the data 
by setting the torsional spring stiffness to K^ = 300 pN- 
fj,ui. To compare Kh to K^, we consider K^/lc, where the 
mean cross-link spacing is given by /c ~ 2 /nm |3]. Thus, 
we find that K^ ~ 20Ki,/lc- One possible reason for a 
larger K^ could be an increase in fiber radius near the 
cross-links by a factor of 2, since bending stiffness scales 
by r^. 

We assume that in the undeformed state of the net- 



work, there are no internal stresses. Thus the fibers have 
an innate curvature and the cross-links have an equilib- 
rium angle equal to that in their initial configuration. 
We ignore entropic contributions to the fiber mechanics. 
While the geometric persistence length of these fibers has 
been measured to be 20 /zm [3], the thermal persistence 
length is much longer Ip = K^/kT ~ 1 cm. Furthermore, 
in the case that the strain stiffening is dominated by ther- 
mal compliance, one would expect to see a decrease in the 
yield strain with increasing concentration |21j . Collagen 
gels, however, have been shown to have a constant yield 
strain of about 60% for a wide range of concentrations 
|22j . Thus the total energy in the network for a given 
configuration is given below. 
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Here Na is the number of elements of type a £ {s, fe, x}, 
which denotes stretching, bending, and cross-link, Li is 
the length of stretching element, Q\ and B\. are the bend- 
ing and cross-link angles respectively, and A indicates the 
difference between the deformed and undeformed state. 

To calculate the stress-strain relationship of our model 
network, we perform a series of 18 incremental strain 
steps by imposing a small deformation on one face, F\ 
while holding the opposite face, Fq fixed. We impose 
two types of deformations: tension (Figure llli) and shear 
(Figure life). In a tensile deformation, we allow the ver- 
tices on Fq and F\ to move freely in directions perpen- 
dicular to the imposed strain to allow for perpendicular 
contraction that is seen to occur in these experiments 
P^ [53] ■ In experiments of this type, the distance be- 
tween Fq and F\ is on the order of centimeters and the 
simulated network represents a small region near the cen- 
ter of a sample. In shear, we do not allow the boundary 
nodes on Fg and Fi to move freely. We compare the shear 
results to cone-plate rheometer experiments, where the 
shear faces are bound to the rheometer. Here, the dis- 
tance between Fq and F\ is 109 /im, and the simulated 
network is one fourth the length of the experimental sam- 
ple between the boundaries. In both deformations, all 
other nodes, including those on the four remaining faces 
of the network, are free to move. 

The minimum energy state of the network at an im- 
posed strain U{eij) is found using a conjugate gradient 
method developed by Hager and Zhang [2^. The stress 
required to hold the network in its current configuration 
is given by (Jij{eij) = {dU/deij)/A, where A denotes the 
area of Fq. The results shown in Figures [2] and [3] are 
averaged over all four extracted networks and over all 6 
principle shear deformations in the sheared network and 
all 3 principle tensile directions in the stretched network. 

The results from our shear cell experiments are given 
in Figures [2] and [3^. In addition, in Figure [HJd, we also 
present the previously reported tensile modulus of large 
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FIG. 2: Stress-strain sweeps of the gel, from 0.5% to 100% 
strain. In the legend, E denotes experiment, W denotes 
the WLC model, S denotes the spring model where we set 
Kx = Kb — 0, and the number denotes the collagen density, 
a) Unsealed results. Note the good agreement between model 
and data at small and large strains, b) When the curves are 
scaled by p^'®** relatively good data collapse is achieved. We 
denote lines of slope m = 1 and m = 3 to guide the eye. At 
large strains, scaling breaks down and the low density curves 
overtake the high density curves because at large strains, stiff- 
ness scales linearly with density. Thus this rescaling serves 
mainly as a visualization tool and does not represent a true 
data collapse. 



samples that are centimeters in length [22^ . Figure l2k 
shows a strain sweep from 0.5% to 100% shear strain in 
both the model and cone-plate rheometer experiments. 
At small strains, the stress-strain response is linear, as 
expected, and at larger shear strains, the stress-strain 
response appears cubic. In Figure [3k, we show that the 
small strain modulus scales by (T12 ~ P^'^^ei2, where p 
is the collagen density. At this time, it is not possible 
to verify the power law scaling in the model since only 
densities of 0.5, 1.2, and 1.4 were observed. The fluores- 
cent labeling of the network changes the polymerization 
properties of the network, causing it to clump at higher 
densities. We use this scaling relationship to collapse the 
curves in Figure [2]d. The close agreement between model 
and experiment indicates that strain stiffening due to the 
geometric rearrangement of the collagen fibers is enough 
to explain the strain stiffening seen in experiments. 

At large strains (~ 50% in tension, ^ 200% in shear), 
the stress-strain curve of the model becomes linear again, 
though with a much steeper slope. In FigurelSb, we com- 
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FIG. 3: a) The small strain shear modulus is compared to 
the cone-plate rheometer experiments and a scaling law of 
G' ~ p^'®* is observed for the experiment, b) The large strain 
tensile modulus from the model and from the experiments of 
Roeder et al. [22]. Results differ by a factor of 2.5, which is 
reasonable since the two experimental protocols were different 



pare the large strain tensile behavior of the model to the 
experiments of Roeder et al. [22]. While our model un- 
derestimates their experimental measurement by a factor 
of 2.5, we find this to be reasonable since the two exper- 
iments used different collagen protocols. In particular, 
different buffers were used. In Figure |2] we also explore 
the case where K^ ^ Kb = 0, such that we have only a 
network of springs connected at freely rotating pin joints. 
At low strains, the network can be deformed without ex- 
erting any stress, but at strains higher than 25%, we see 
that this simplification adequately describes the gel. 

A topic of investigation explored by many is the valid- 
ity of the assumption that these networks deform aflinely 
[TTl E5] . For brevity, we state only that the deformations 
are highly nonafline in these simulations, as evidenced 
by a visual inspection of Figure IT}! and Wp where it is 
obvious that many fibers leave the volume defined by an 
affine deformation. 

In summary, we have presented a microstructural 
model of a 3d biopolymer gel using a network geome- 
try that is based on the true network architecture. It 
differs from previous work in that we use realistic net- 
work architectures that have been extracted using the 
FIRE algorithm. We specifically focus on the mechanics 
of coUagen-I networks, but emphasize that this modeling 
approach is generalizable to other biopolymer networks. 
The model has three parameters: {E,r,Kx}. The fiber 
radius and tensile modulus can be measured experimen- 
tally and the model uses realistic parameters. The cross- 
link torsional spring constant must be fit to the data and 
we used K^ ~ 20Ki,/lc- Fitting this single parameter 
gives the right strain stiffening behavior for networks at 
three different densities at strains that vary from 0.5% to 
50%. This result lends support to the hypothesis put for- 
ward by Onck et al. [llj that strain stiffening in polymer 
networks, and particularly coUagen-I gels, is governed by 
rearrangement of the gel. 

Another finding of the model is that at strains greater 
than 25%, the stiffness of the gel is governed almost en- 



tirely by stretching of the fibers. This result is relevant 
for collagen because cells embedded in these gels are seen 
to produce deformations of this order of magnitude [26] . 
In modeling large systems of this type where the strains 
are large, it may be sufficient to treat each fiber as a 
spring rather than a WLC in order to reduce the compu- 
tation time. This work also demonstrates that an under- 
standing of the cross-link mechanics in these systems is 
critical to understanding their mechanical properties, as 
has been seen previously f27; . In much of the theoretical 
work that has been done on random stick networks, the 
cross-links are treated either as freely rotating pin joints, 
or welded joints of fixed angle jTU [2TJ [25]. While these 
are sensible simplifying assumptions in developing a the- 
ory, they are not adequate for describing actual networks. 
We note that this model has been designed to capture 
the short time scale (< 1 hr) behavior of the network, 
where it behaves as an elastic solid. Such behavior re- 
quires that the cross- links be relatively fixed [28]. This 
simplified model provides a starting point in the devel- 
opment of a more complete model of collagen gel. Ul- 
timately, a more sophisticated approach, such as that 
taken by Rodney et al. [29 will be necessary to capture 
the full dynamic behavior of the gel, where cross-links 
are allowed to slip and break. 
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